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ABSTRACT 

The thesis covers the mechanical, thermodynamics, and aerodynamics 
of designing a free vortex axial flow turbine with the aid of a digital 
computer, It uses the Improved Ai nley-Mat hi eson Performance Estimation 
Method to determine the aerodynamic pressure losses so off-design and 
design point calculations can be conducted to produce performance curves. 
It also gives physical and thermodynamic properties of the turbine. 
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The purpose of this thesis is to develop a computer program 
which can be used in the Preliminary Design of a Free Vortex Axial 
Flow Turbine. The design plan consists of three main parts. The 
first will be cycle calculations to determine stage characteris- 
tics at each station such as angles, density, temperature, flow areas, 
etc. The properties from this phase will be used in parts 2 and 3. 

The second part is to determine mechanical characteristics such as 
size and weights of turbine parts and stress and vibrational char- 
acteristics. The final phase is to determine turbine characteris- 
tics at design and off design parameters such as speed, pressure ratio, 
etc. This last phase is optional since compressible flow tables and 
gas properties are required and many turbine designers may be more in- 
terested in the size and weights than off design characteristics in 
the preliminary design phase. These three parts are covered in cha- 
pters one through three, chapter four will cover the operation of the 
computer program developed in this text. Appendix one is a list of 
symbols used in the text, appendix two is a user f s guide to the com- 
puter program and has an example turbine design which is compared 
with a design proposed by reference 20. Appendix three is the com- 
puter program. 

The computer program is set up so it has many options of in- 
put variables so that as the designer becomes more knowledgeable 
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in his design study he can specify more characteristics and have 
more control over the design output. There are checks in the program 
to limit poor design properties such as excessive swirl in output, 
high tip to hub ratios, negative reactions, etc. The computer pro- 
gram can be used to design a constant hub, mean, or tip turbine. The 
design does not consider the effects of cooling in turbines, since 
time was not available to incorporate this into the computer program. 
This limitation mostly affects aircraft gas turbine and aircraft de- 
rivative marine gas turbines since these engines operate at high tur- 
inlet temperatures and cooling is required to allow for a high mean 
time between failures. 

There has been no effort made to prove the relations used in the 
design calculations. These can be found in basic physics, thermody- 
namics, and fluid mechanics books. 

The notation as to signs of angles in the velocity diagrams will 
be the same as used in chapter 4 of Gas Turbine Handbook which is re- 
ference (3). An illustration is enclosed in chapter l of the sign 
convention and notation used in this text. A deviation from this no- 
tation will be used in section (3.2). 



Chapter 1 
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Cycle Calculations 



1.1 The basic design of an axial flow turbine starts with the 



cycle calculations to determine the thermodynamic properties thro- 
ughout the turbine. The basis of these calculations are from known 
inlet conditions and assumed design parameters. The computer pro- 
gram is set up so that a minimum number of these parameters must be 
specified and leaves the option to the designer to more fully spec- 
ify the design parameters, if he so desires. This chapter intends 
to show how these basic calculations are carried out and the methods 
employed to determine the physical properties such as flow area from 
these design parameters. Figures (1-1), (1-2), and (1-3) show the 
sign convention and notation that will be used in this chapter. The 

following are some of the design parameters needed for the cycle cal- 

o 

culations; inlet temperature ( R) Toi, turbine inlet pressure (psi) 
Poi, omega (rpm)/2_, specific heat (btu/lbm °R) Cp, mass flow 



( 1 bm../sec)m, ratio of specific heat gamma 6 . tip radius (ft)r t ,tip to 



hub ratio (rt/r^ ) , 




etc. 



Appendix two (2) will specify which parameters must be known and 



which ones are optional. 



1.2 Determination of overall turbine parameters from inlet 



conditions and assumed design parameters. 




( 1 . 2 - 1 ) 



( 1 . 2 - 2 ) 



Toe = Toi - ^ hot (°R) 
Cp 
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Fig. l-l Axial Flow Turbine Stage Fig. 1-2 Velocity Diagram 




Fig. 1-3 T-S diagram for a Reaction si.v-e 
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Polytropic Efficiency 

ln(Toe/Toi ) 



% 



In 



( Toe/ Toi- 1 ) + i 

7Ct J 

r 

Poe = Poi ( Toe / To i ) ' f\p 1^-0 



(1.2-3) 



(1.2-4) 



(b) If the pressure ratio were specified insteaa of ootput power 
the following calculations will be performed. 

Polytropic Efficiency ^ 

0 U. UO -ftO-CPoe/Pon 4 ^) 

/L r sy 

ln((Poe/ Poi )-^>;.-j“— _j 



y 

Toe = Toi (Poe/Poi) 9^ p (V- 1) (1.2-6) 

/\h° = Cp( Toi - Toe) (1.2-7) 

Power out = . (1.2-3) 

550 

(c) 

C ( rad/ sec) (1.2-9) 

60 

r h = r t /(r t /r h ) (1.2-10) 

r m - (r h + r t )/2 (1.2-11) 

Uh “ r^&/( ft/sec) (1.2-12) 



Euler’s Equation 



Ahc 



f 



stage 



UhVx (taiflfe - tan(X^) 
g° J 



^oJAho st:a ^ { 

u h 2 
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(1.2-13) 

(1.2-14) 



The value of l^is taken at the hub of the last stage because this 
is the most critical in the design of a turbine in reference to load- 
ing. Normally the value of^is from (1.5 - 2.9). The lower value is 
for lightly loaded stages while the higher is for highly loaded stages, 
To have no exit swirl must be less than or equal to 2.0. 



l/)^ 2.' 



( no exi t swi r 1 ) 



for^~ 0 



| 2 



where 



(1.2-15) 




h E. - ^3 
ho^,- hob 



ratio of change in static enthalpy 
across the rotor to change in 
stagnation enthalpy across the stage. 



With the value of 




speci f ied , 



the number of stages required to 



develop the power out desired can be determined. 



A ho 



stage 




( 1 . 2 - 16 ) 



Determining number of stages 



N 



integer 



Aho t 

Aho sta ^ e 



If N were specified instead of 




the change in enthalpy per 



stage would be determined. 



A h ° 



stage 



Z)ho t 



N 



If 



This program uses a 

'^2.0 and(Xh " 0 

equation (1*2-17). 
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constant change in enthalpy for each stage, 
the radius of the huh can be determined from 



r he 



goJ Aho s tage 

U> 2 



(1.2-17) 



Another parameter which will be used is ^P^Vx/U^ (flow coeffic- 
ient). Where Vx is the axial velocity, which is constant in a free 
vortex turbine. This parameter ranges from 0.4 - 1.3, the higher 
values are used for aircraft turbines which use high axial velocity 
for propulsion. In marine and industrial turbines a low value is 
desired because the energy of the exhaust gases are not recoverable 
in this form, unless it is used to drive another turbine. 

Another parameter that may be specified is the critical Mach 
number at the exit. 



Merit 3 



Vx 



goRToe 



where R is the gas constant 
R = CpO^l/tfj 
The Vx is computed exactly 




)) goRToe 



0 * 



(1.2-18) 



(1.2-19) 



Vx - M 



( 1 . 2 - 1 8 a ) 
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If the exit Mach number is specified the turbine exit area is 
varied until the required area is obtained to give the specified ax- 
ial velocity Vx. If this number is not specified but the exit area 
is specified the value of Vx is varied until the required area is 
achieved. 

From the continuity equation for gases 
Vx =» ™ 



where 



2/>/^r)rdr 



f (r> 



h 

P(r) 

RT(r) 



T(r) = To - 



V(r)‘ 



2CpgoJ 



Tf" 

P(r) = Po/ ( To/ T( r ) )T-T 



( 1 . 2 - 20 ) 



( 1 . 2 - 21 ) 



( 1 . 2 - 22 ) 



( 1 . 2 - 23 ) 



V(r) - 



Vx 



cos®^(r)) 

where r varies between the hub and tip. If the exit angle is not 
equal to zero 



( 1 . 2 - 24 ) 



-1 



O -$he - fn >/ jrf 



O^Viubg^ tan 
and for free vortex 

rVe^ r Vx( tan(y ) 5 constant 

then 






ton 



-1 



r he tan ( (Vhe ) 



( 1 . 2 - 25 ) 



( 1 . 2 - 26 ) 



( 1 . 2 - 27 ) 



where re is taken at hub mean, and tip. This calculation is carried 
out in function program ANGLE. Thus the exit angle has been deter- 
mined at the hub, mean, and tip radii. 
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The pressure, density, and temperature calculations are per- 
formed in subroutine TPl . The integral calculation is performed in 
subroutine DENR which uses Simpson’s Rule to compute the axial vel- 
ocity Vx. With the absolute angles 0C( r e ) determined the relative- 
angles (r e ; can then be determined by equation (1.2-28).. This is 
carried out in function program ANGLES. 



S( r e ) = tan 
From equation (1.2-15a) 

■ r j * 1 . -4ho s tage^°^ 






tan(6<l(r e )) - 



U(r e ) 



Vx 



2(CUr e ) 



Vxtan(04) 

0)r 



(1.2-28) 



(1.2-29) 



which is computed in function program REACH. With the exit cond- 
itions specified, the properties at station 1 and 2 can then be 
determined for each stage. 

The following properties are now known at the exit of the turbine; 
Vx, Toe, Poe, J ? (r g ), P(r g ), T ( r e ), ^( r e ), f£) ( r e ) , (5V r e ) » and r e * where 



1.3 



r, , , r , r . 
hub’ mean tip 

Determination of properties at station 1 for stage n where 



n — 1 . . . . N. 

To, = To 3n + (1.3-1) 

CP 

r 

Poln - Po3 n /(To3n/Tol n ) ^fT) (1.3-2) 

If this is the 1st stageO^/n = 0.0 where n = 1 

If not (X/n( r c ) = £)<^ n (r c ) where r c is the radius which remains 
constant through out the turbine, ie(constant hub, tip, or mean). 






tan 



tan<(y (r (J )) - U( 



Vx 



T/n(i c ) *» To/n 



v lf£ r c) 

2goJCp 



V in('-c) 



Vx 



cos (Oh n( r„ ) ) 



Y 



Pjn(r c ) = Po/ n/( To/ n/T/n) Y") 
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(1.3-3) 



(1.3-4) 



(1.3-5) 



(1.3-6) 



vj 



«n(r c ) - 



(1.3-7) 



RI>n(r c ) 

which is identical to equations in section for exit of the turbine. 



A trial blade height (h) at station In will be established. 



h = 



m 



(1.3-8) 



2?Tvxr c ^(r c ) 

Let delta h = h/lO.'-'Now depending on whether the turbine is constant 
hub, mean, or tip, radius are established at the other radii. 
Constant hub; 

r ,n(hub)= r c (1.3-9) 



r /n (mean) « r c + h/2 (1.3-10) 

r jn (tip)» r c + h (1.3-11) 

Similarly for constant mean or constant tip. 

can then be determined at all radii for station 1 and temp- 
erature; pressure, and density can be calculated at all radii in the 
same manner .as was carried out at the exit of. the turbine. Then 



Vx prime cart be calculated in the same manner as Vx was calculated 
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Vx prime 



ft 



(1.3-12) 



2(7'] b?r)rdr 
•/ r h a 

When Vx prime = Vx the proper radii has been calculated at station 
l. The method employed continuously decreases delta h until conver- 
gence is accomplished. This is the same method which was used at the 
exit of the turbine and will be used in station 2 calculations. The 
properties of the gas is now determined as per exit of the turbine for 
station 1. If tnis isn't the first stage the properties at station 1 
are the same as station 3 of the preceeding stage. 

1.4 Determination of properties at station 2 for stage n where n = 

1 . . . . N 

To2n = Toin (1.4-1) 



(XtvS- r c ) = tan’ WjAhOs * tan(<^n(r c )) 
^ Vx cOr c 



(1.4-2) 



V(r ) =• Vx/cos(^ 2 n(r )) 



(1.4-3) 



T 2 .n = Ton.n - V(r c )/2goJCp 






(1.4-4) 



Pzn p ( n( T^n/ Tjn) ->lp()f-l ) (1.4-3) 

which is an approximation from figure (1-3), where the actual pressure 
is determined from pressure 1,^*$ data as will be covered it Chapter 3 . 



POW ~ P2.a( ro2n/ T^n/cP" ^ 



(K4-6) 






P2Q(r ) 

HTzn( r ) 



(1.4-7) 








\ 
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These calculations are accomplished in subroutine TPD2 to 
determine temperature, pressure, and density at station 2. The calcu- 
lations are carried out in the same manner as station 1 calculations f 
to determine the gas properties and radii at station 2. The calculations 
are then carried out at station 1 of the next stage(ie; n = n-1 until 
properties at all stations and positions have been computed for all 
stages, ) 

1.5 Having solved for gas properties at all stations and positions 
and having flow areas and gas angles the relative stagnation temperature, 
and the relative Mach numbers for stators and rocors can be calculated. 



To2rel~ T + " ( i • 5 - 1 ) 

2goCpJ 









-> 


2 






M nozzle inlet 


= 


goJn/Tjn) - I 


>r-i 






(1.5-2) 






v_ 














r 




h 






M nozzle exi t 


1 


poln/Tzn) - 1 








(1.5-3) 






L 














r 






\h 




Mrotor 

i nlet 


=2 


^To2ji rel /T2.n) 


-3 


y '*j 


\ 


(1.5-4) 






r 






h 




M 

rotor eX i t 


as 


((To3n rel /T3n) 

>s_ 


-3 


2 

J 




(1.5-3) 



These calculations are carried out at each stage. The following 
parameters have been determined at all stations of the turbine; temp- 
erature, pressure, density, flow area, absolute and relative angles, 
Reaction, relative Mach numbers, and the axial velocity Vx. These 
properties will be used to determine the mechanical aspects of the 
Turbine design anc! be used in performance estimation. 



Chapter 2 
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Mechanical Design of an AxiaL Flow Turbine 

2.1 In this chapter the mechanics of the program and assumptions 
made to design nozzle and rotor blades, and the turbine disc will be 
discussed. 

2.2 Having the gas properties and flow angles throughout the tur- 
bine the blade angles can be calculated for the turbine. The criter- 
ia used in establishing the blade angles is the loss data provided by 
figure (2-2) and (2-3) based upon Ainley Mathieson method of perform- 
ance estimation. As can be seen for reaction blades in figure (2-2) 
the incidence angles can vary from -15 to + 15 without much of an in- 
crease in profile loss Tp. In actual construction this large varia- 
tion can allow for a decrease in twist in the blade from hub to tip 
normally required in free vortex blades. For the computer program, 

selection of blade angle, at the inlet of the blade is for zero incid- 
ence at design point. 

Figure (2-3) shows how blade exit angle varies with gas angle 
there is also a correction for Mach number. For Mach numbers less 
than 0.5 the values given by figure (2-3) will be used, for values 
greater than Mach 1 the blade angles will be equal to gas angles 
while for values between Mach 0. 5-1.0 a linear variation will be as- 
sumed as is shown in figure (2-4). The calculation of these blade 
angles are carried out in subroutine 31ade, which also will be used 
to determine gas angles when blade angles are known for off-design 
calculations in chapter 3. There is an additional correction for blades 
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Figure 2-1. Conventional Blade 
Prof i ie 



Figure 2-2. Variation in profile 
loss with incidence angie(Fro:n 
D.G.Ainley and G. C. R. Ma thieson( 1 ) ) 




Figure 2-3. Relation between 
gas and blade angles (From 
D.G.Ainley and G. C. R. Mathieson( 2 ) ) 




0.0 0.23 0. 5 0. 73 1.0 1.23 

Outlet Mach Number 

Figure 2-4. Angie variation 
with mach number ( From D.G.Ainley 
and G.C.R.Mathieson(2)) 
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with a curved back trailing edge in reference (2) but this will not 
be applied in the computer program. 

2.3 The computer program can now determine the ideal pitch to chord 

ratio P/C, blade chord (C), aspect ratio H/C, bending and centrifugal 
stresses, number of blades, 1st bending frequency, and weight of blades. 
As will be shown in this section all these parameters are interrelated 
and they all help determine each other. 

2.3a Centrifugal Stress Calculations 



dr 

Ab 




4 

r 



r 



t 



r h 



cO 




d F = dmr CO^* = dmr (2TMX/60) 



Figure 2-5. Rotating Turbine blade 
2 



(2.3-1) 



A(r) * Ab f(r) 



dm 



Ab f (r 



)dr t|r 



dF = ( 277#/ 60)— A(r)rdr 
If A(r) * constant 
Fc 



(27K1/60) 



r ti P 



r hub 



<^ C 



— - (2 7ra/60) 2 -w£- 

Ab go 



(2.3-2) 

(2.3-3) 



(2.3-4) 



2 

r tip 



2 

r hub 
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1.0 



Taper factor 

ff 0.8 



0.6 



0.4 



0.2 



V 

Linear Taper iP' 

7 . // I 









. /' 



. , Conioal lapjer 

V 



/ S' 

■'y/ 



/ 



'/ i 



i 

f- 



-L- 



0.0 0.2 0.4 0.6 0.8 1.0 

Area Ratio - Tip Profile Area 
3ase Profile Area 

Figure 2-0. Variation of taper factor with oladc 
profile area ratio (From Emmert, Current Design 
Pract i ces_£pr Qas Turbine Elements . (9 ) ) 



(2.3-5) 



c - 4.51 Tf 0 Ab LQr-\ 

J Viooo/ 

where; *» average tensile stress(psi) 

Tf - taper factor(from figure (2-6)) 

* specific mass of the blade 
A - annular flow area of blade ring (sq. in.) 

XL 

- shaft speed (rpm) 

For the computer program Tf will be based on linear variation with 
tip to hub area ratio (areara) as an input, normally this will be 



be tween 0.2 5 - 0. 333 



If a linear taper is assumed equation (2.3-1) can be integrated 
directly assuming; 

A(r) « Ab(l -0(r) (2.3-6) 
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r c * 



^AbJo‘-Xr)rdr 
rhub 

2 

2ira\ /b 



(2.3-7) 



L )/ b Ab (rtip - r^ub ) _0^ , 3 3 ) ! 

60fi%o Ab L 2 3 ' r hub r tip 'J 



) _ I — ( r ^jp ~ l , 3 3 71 

©c " Ab = ( 60 / go|_ 2 + 3 ' r hub t i p* IJ 

and can be evaluated 

A t ip - OC area ratio (Ab) 



(2.3-8) 



0^ 



1 - area ratio 
r tip 



(2.3-9) 



2,3b Determination of ideal pitch to chord ratio. This value for 
P/C ratio is taken from figure (2-7) which is a correlation of Ainley's 
optimum spacing data. For stators the inlet gas angle to enter figure 
(2-7) is - and the relative efflux angle is For rotors ^ n is 

the inlet gas angle while -^^is the relative efflux angle. This cal- 
culation is done in subroutine STOCRA and the input angles are taken 
at the mean radius. 

2.3c Determination of non-dimensional section modulus(SM) and non- 
dimensional area at the base of the blade (Area/C^). An input to get 
these values is the turning angle at the base of the blade. ie(hub 
for rotor and tip for stator). Figure (2-3) taken from referenced) 
illustrates how the non dimensional base area is obtained. The non- 



dimensional section modulus is obtained from figure* (2-9) which is 
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Figure 2-7. Optimum spacing of turbine 
blades. (From a correlation of Ainley^ 
loss data (2)) 



an unpublished relationship provided by Ainley for use in reference (6). 
A modification was made to this relationship based upon calculations 
performed by Professor A. D* Carmichael to make the section modulus 
larger to conform to modern turbine blades. The SM was multiplied by 
2.35 for use in the computer program. 

2.3d Determining bending stresses on turbine blades. As can be seen 
in figure (2-10) and equation (2.3-10) a detailed calculation would be 
needed to determine the gas bending stresses on a turbine blade but 
because angle ^ is small and Mw is by far the greater bending moment an 
approximation proposed in reference (6) is used in the preliminary 
design. 






max 



ft Vx ( t a n(X 2 - tanc<J 

- — - 

2 go n Z 



( rotor) 



(2.3-11) 



v^here the angles are taken at mean radius and the following are the 
definition of the parameters; 



0.24 




20 40 60 80 100 120 



Blade Turning Angle 

Figure 2-8. Approximate Areas of Turbine 3lades vs 
Turning angle. (From paper by R.E.Dundas in Gas 
Turbine Handbook(7)) 




Figure 2-9. Section Modulus as a function of blade 
turning angle. (From unpublished paper by Ainley 
for use in reference (6)) 
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X 




C$ab " — (Ma cos f - Mw sind>) 
* lyy 

+X_( Mw cos <£ ♦ Ma s i n <£ ) 

I XX x x 

(2.3-10) 



y 



Figure 2-10. Sending moments acting 
on a turbine blade. 



L - SM C Section Modulus 

SM * non dimensional section modulus 

n - number of blades 

* 2 *n~ r m where r m * mean 

Pi tch 



radius 



(dgb) mo .»“ 



m Vx (tanflC?- tanOCg) h 
2 go 2^ SM C 3 



(2.3-lla) 



ft Vx ( tan 0C2 - tan0(3) h P 
" ATIlor m SM C2 C 

At this point the blade chord has not been determined so the bending 
stress cannot be calculated, but equation (2.3-lla) wiLl be coupled 
with a natural frequency equation and a cyclic loading equation to 



calculate the blade chord. 
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Figure 2-11 Effect of taper on blade natural 

frequency ( Frorr Emmert, Current Design Practices 

for Gas Turbine Power Elements (9)) 




Blade Height Ratio 



31ade Height 

Blade Ring Mean Diameter 



Figure 2-12. Variation of rotation coefficient with 
blade height ratio. (From Emmert, Current Design 
Practices for Gas Turbine Power Elements . (9 ) ) 



2.3C 



Determination oC natural bending frequency of blades from 23 



reference (9). The following expression is used in the computer 
program. 



f s = 11.0 C f 



k 

"h 7 



(2. 3-1 2a ) 



>ljf 



f s is standing frequency (cps) 



Cf frequency correction factor taken from figure (2-11) 
k minimum radius of gyration of base prof ile( Iuiu) 
h blade height (in) 

E modulus of elasticity (psi) (^29 - 30]x 10^ psi 
J>b specific weight of blade material (lb/in-^) 

f 



-s - 11*0 % 
h“ 



S 



IE 



>*•<*& 



ZtE 



■S 



assuming I t 



A cross sectional area of base of blade (in^) 
t maximum thickness of blade (in) 

- (t/C) C 
A = (A/C 2 )C 2 



n.o Si 

h 2 



n.o Si 

n 



x 



SM C3(t/C)E C 



(A/C 2 ) 



SM(t/C)E 



(A/C 2 ) _ 



chord 



(2*3-1 2b) 



( 2 • 3- 1 2c) 



The running frequency f given by equation (2*3-13) is a function 



of the natural frequency and the speed of rotation of the blade. The 
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factor d is given by figure (2-12). A plot of this frequency as a 
function of operating speed is illustrated in figure (2-13) which 
is called a Cambell Diagram. Normally for good design practice the 
running frequency should avoid the 1st four harMBrtfcs erf the operating 
speed in its operating ranga. Ear the design of the blade chord :in 
the computer program the ief feet o± the rotational speed as indicated 
in equation (2.3-13) will be neglected but a more stringent criteria 
uding the' harmonics of. the: operating speed to estimate the ivibxlation- 
di stress in the bLade will be used. 



f 



r 



Jf2 + 3(-^-/60) 2 



( 2 % 3 - 1 3 ) 
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Frequency 

(cps) 
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Figure 2-13* Cambell diagram, plot of frequency 
of vibrations as a function of operating speed. 



Detertnl nation of* turbirte\ blade dhAlracteriatlca bA&ed on 
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2* 3£ 



cyclic loading requirements by use of a Goodman Diagram of steady 
and vibrational stresses. In reference (7) a use of the Goodman 
diagram is used to determine the chord of a compressor blade. This 
can also be applied to turbine blades. A Goodman diagram is a plot 
of vibrational stresses as a function of steady stresses. The vibra- 
tional stress is taken at 10^ cycles at the operating temperature. 




Figure 2-14. Goodman Diagram of Vibrational 
stresses as a function of steady stresses. 

This plot can be approximated as a linear function. 




In reference (7) the vibrational stress is proposed as a function 
of the bending stress as shown in equation (2.3-14) 



vib = 1.3 Amplification factor(3>§b 



(2.3-14) 



31 




Figure 2-16. Amplification Factor. (From 
Paper of Trumpler & Owens (19)) 



where Amplification factor is shown in figure (2-16) as a function 
of the harmonic of operating speed. The 1.3 is a stress concentration 
factor specified in reference (7). From personal conversation with 
Professor A* Douglas Carmichael at M. I.T. the amplification factor 
can be approximated by 41/n where n is the harmonic of the operating 
speed. 



n 





60 



(2.3-15) 



From figure (2-15) 

^vib - Cl - £i^steady (2.3-16) 

where ^ steady * ^.centrifugal + ^ gas bending 
and Cl is. the upper limit for vibrational stresses and C2 is the 
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upper Limit for steady stresses. 




1.3°/6 0 i CJ_ 
f C2 




but 



a 

C2 



^ 1.3 (-0-/60)41 



'O ( 1 . 3x41 x P-/60) ri/'s \ 

SO Q ?;> f ^ C1 " ^ Cf ) 



(2 



From equation (2.3-lla), (2.3-12c), and (2.3-17) 



mVx( tarfrC? - canfrff) /p\h( 1 3x4 lx(l^/ (?Q ) _ ^ 

4ffeo r_ SM ? VC lL.OxCr 



SM(t/C)£ 



_(A/C 2 )j 



C 2 



(2 



Stress 3 = 



:nVx(tarC(2 - tanX3)h /pj 
AiTgo r m SM \C/ 



(2 



Freq = Ll.OxCf 
h 2 



SM(t /C) 

_(A/C 7 )j!b 



(2 



^^ Al CiVo O) (Stress3) _ cl(1 .6c£) 



"T 

C J Freq 



C2 



(2 



C4 = 1. 3x4 L (£1/60) Stress 3/ Freq 



(2 





f-Scf 
C2 1 



.3-13) 

.3-19) 

.3-20) 
.3- 13a) 
.3-21 ) 



Chord 



1/3 



(2.3-22) 
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Now that the chord has been determined the bending Jtress and natural 
bending frequency can be determined from equation (2.3-lla? sund (2.3- 
12c) respectively. Having solved for the chord at the base of the blade 
the chord can be determined at mean and tip assuming a linear taper. 

(2.3-23) 

Assuming constant (A/C^) 

(2.3-24) 



, , Ab( 1 + Areara) 

A Dm = — 



Cm = Abm ^ * 

jaTF) 

At the opposite end from the base 

r _ L\rearati o( Ab) 

L t \ : — ~2 



L (A/C z ) _ 



(2.3-23) 



where C c = C hub 



( for stator) 
C t = C^. 5 n ( for rotor) 



C t ip 

Now that the chord has been determined at the midchord the aspect ratio 
and pitch can be determined by equation (2.3-26) and (2.3-27). 

Aspect ratio (h/C) « h/Cm (2.3-26) 



Pitch m - (P/C) Cm (2.3-27) 

With pitch determined the number of blades can be calculated, 

2Tprm 



Number of blades ( N) =* Integer 



(2.3-23) 



Pi tch 

The computer program uses the criteria having even number of blades 
for the stator and prime numbers for the rotor to prevent blade inter- 
action of resonant frequencies. The determination of whether a number 
is prime is done in subroutine TCSTP. Thus havin' solved for the 



number of blades the pitch can be recalculated and the weight 



of the blades. 
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Blades are integrated from the base to end to get 
volume and is multiply by its specfic weight to get the total weight. 
In the computer program this specific weight is used for both the 
blades and disc. Also it is assumed they are solid blades with no 
cooling passages. 



Weight of blades 




(2. 3-29) 




where Otis given by equation (2.3-9) 



r h > 





In this section pitch was determined by a chart of ideal pitch to 
reduce losses. In actual turbine design the loss data would be 
coupled with stress analysis of the turbine blade attachment to 
determine allowable pitch. The tyon nost commonly used is the fir- 
tree method which is illustrated in figure (2-17). This method of 
selecting turbine blade chord length is just one of many which are 




Pi 

at 



Figure 2-17. Firtree olnde 
attachments . 
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2.4 Turbine Disc Design 



The method for design of a disc in this section will be based upon 
a constant average stress analysis. The stress relations are those 
proposed in reference (15) and (18). This method will ignore stress 
concentrations ana thermal stresses. In actual disc design these cannot 
be ignored since thermal stresses may be a major part of the disc 
stresses and stress concentrations may cause the disk to fail even th- 
ough the average stress is well below the ultimate tensile stress of 
the material (UTS). For calculations in the computer program the aver- 
age tensile stress used is that shown in equation (2.4-1) from reference 



The value of this ratio is normally between 1.2 - 1.3. The value of 
ultimate tensile stress (UTS) is taken at the design operating tem- 
perature, higher values producing smaller disc dimensions. This 
value is an input to the computer program. The computer program has 
minimum values for some of the dimensions illustrated in figure (2-13) 
to prevent impossibilities to manufacture and to allow adequate thick- 
ness for attachment of blades 

The following equations are based on figure (2-13) taken from reference 



(7) 




(2.4-1) 




Ratio of burst speed to design speed (2.4-2) 



( 1 ." )• 




(2.4-3) 




Wo To 
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c<b - 



1 - 



r ii 



20{r t 



2go 



1 - 



! rh> 



(2.4-4) 



for a linear tapered blade. 

For the computer program Wo will be assumed to equal .65 chord 
and To will equal half of Wo unless otherwise specified. With these 
parameters known Zo can be solved for by equation (2.4-5). With Zo 



Zo 



/ Wo To 



£o 

LVt, 



\ 2 + ^ 

' 217 Wo To 



i2_ 2 



■P ufrt 

T* 

2dgo 



(2.4-5) 



known Z(r) can be solved for. 

2 



■ f 

Z(r) = Zo e' J S 2go 



£a 
rt/' 



,V. 



(r/J 



(2.4-6) 



with a minimum of 0.3 inches. The value of r^ is assumed to equal 
2 r s | ia £ t where r s haft is an in P ut to t ^ ie computer program. This 



defines Ti so Wi nn be solve for in equation (2.4-7). 



Wi = 



Zir j 

Ti 



1 - M 2 fU) 2 r t 2 

W 



(2.4-7) 



with a minimum of 1.2 Wo. 

Having determined the dimensions of the disc the weight of the 



disc can be determined by equation (2.4-8) 






Disc weight 



- j/i 

o r 



A( r )drdU 



(2.4-3) 



A(r )dr 



i c tL c r ° r Vh 

2 TT Wi J rdr + ) z(r)rdr + Wo I rdr 






= 2l>j>f^( r i 2 - r shaft 2 ) + Wo( r h 2 - r Q 2 ) 

^ n ,. x 2~ 2 r „ 2- 2 ,_.x 2 _ .2 2 



+ 2Zor t e 



■PC p'-fo 
2 3-^go 



f -pu^r t 2 / ri \ 2 -raJ 2 r t 2 /r 0 \ 2 1 

|e - e (jT,) J 
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Having determined the dimensions of the disc, the burst speed 
can be computed assuming the disc yields prior to bursting. For this 
analysis the disc is divided in half as illustrated in figure (2-19) 
and the stress is assumed constant over area 2A. The force acting 
on this area is composed of the disc load and the blade load. It 
is assumed that the blade load is evenly distributed on the outer 
rim. In actual stress analysis of the disc this assumption is an 
over simplification of the stress, but for preliminary design 
purposes it will be adequate. A more detailed analysis is covered 
in reference (7). The force due to the disc is now calculated by 
taking the force in the direction perpendicular to the area 2A 
caused by the disc load and summing it over the disc half and adding 
this to the force caused by the blade load. This force is divided 
by the area 2A to get the average tensile stress^ ta. This stress 
is set equal to the ultimate tensile strength of the material of the 
disc at its operating temperature. The value of Nb is then calculated 
by equation (2.4-16) 
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dF - 



dF 



di-c 

+ ^blades 






Figure 2-19 
Turbine Disc 
load distribution 



A 3 Wi Ti + 


r ro 

^Z(r)dr 


For disc load 





+ Wo To 



dm 




V 



(2.4-9) 



(2.4-10) 



dF - (dmro?)sinW 





(2.4-11) 



(2.4-12) 
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_ f- a r 2 ? 



if / , \1 




<o 

•'av'. 



v go L> r -t' 



e^tc(due Co disc) = ■ ^ ^ ' ■ i 



(2.4-13) 



For blade load 



Qj}ta(due Co blades) = ~ <^b given by equac i on( 2 . 4-4 ) 



(2.4-14) 



Average disc stress 



C$ta(Disc + blades) = r dise + ^blades 

2A 



(2.4-15) 



Set c^ta = ULT (ultimate tensile strength of material) 



The burst speed Nb can now be calculated. This value will 
slightly higher than that in dquation (2.4-2) due to approximations and 
arbitrary values selected for some of the disc dimensions. 



Nb = 







d 



Ult 



2 A 



■S 2 f. 



^(r)r“dr + 



r shaf t 



NAbr 
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u/rhV- ^ rL 
. bey ^ 



i 



i 




(2.4-16) 
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2.5 Having determined the blade and disc dimensions the axial 

length of the turbine can be calculated. To do this a criteria for 
spacing between blades must be selected. From reference (6) the 
values of spacing between blades is 0.2 to 0,5 of the axial chord 
length. The longer dimension will tend to reduce the flare which 
would prevent the gas flow to separate from the turbine wall. It 
would also decrease the chance of inducing vibrational stresses. 

The shorter dimension would decrease the size and weight of the 
turbine. For the computer program a value of \ the average axial 
chord at mid-height was selected. The computer program will then 
calculate the axial length and weight of the shaft from station 1 
of the 1st stage to station 3 of the last stage. A value of 13° 
was selected as the stagger angle in the computer program. 




Figure 2-21 Stagger Angle 
of Turbine blades 



Chapter 3 
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Performance Estimat ion 

3.1 The design and off-design characteristics are very important 

in a complete analysis of a gas turbine since most turbines do not 
operate only at their design point. Also cycle calculations without 
actual loss data can not accurately determine the efficiency of a 
turbine, so this analysis can give a more accurate assessment of the 
efficiency that can be achieved in the turbine. 

The Ai nley-Mathieson Method for performance estimation was 
selected over other methods because it seems to have a wider applic- 
ation and can be easily used for off-design performance calculations. 

The modifications proposed in reference (8) have been applied. These 
modifications can reduce the error fromi‘3/o to t 27* for most turbines 
and are based on modern turbines while the data by Ainley and Math- 
ieson were accumulated over 20 years ago on turbines that existed at 
that time. The calculation procedure for off-design performance is 
that proposed by Horlock in reference (12). A description of the 
Improved Ainley Mathieson Performance Estimation Method will be pre- 
sented first followed by a description of off-design calculations which 
use loss coefficients to determine gas properties, isentropic efficiency 
and static efficiency at design and off-design points. 

The computer program used pressure loss coefficient Y * ( Po: - PoC ) 
/ (P 02 .-P 2 . ) for performance analysis. This coefficient is a func- 

tion of aspect ratio (H/C), pitch to chord ratio (P/C), velocity, Mach 
number, trailing edge thickness, thickness to chord ratio (t/C), 
incidence angle, tip clearance, Reynolds number, etc. This coefficient 



A3 



* 



Loss Coefficient 
Y - p °l -P02- 
P 02 -P 2 




Figure 3-1. Analysis of losses in flow 
through a row of turbine blades( From 
D.G.Ainley and 'l. C. R. Mathieson(l ) ) 




Nozzle iilades fa = o> 




Pi tch/Chord( P/C) 

(a) 



Impulse Blades (fit = OC 2 ) 




Figure 3-3. Profile-loss coefficients for conventional 
section blades at zero i ncidence. t/c - 20 per cent; Re 
- 2x10*; M Q^kiFrom D.G.Ainley and G. C. R.Mathieson(2) ) 
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is broken down into three parts, Vp (profile loss), Ys (secondary loss), 
and Yk (clearance loss). Figure (3-1) from reference (1) shows the 
first two, Yp and Ys in which the annulus loss is grouped into second- 
ary losses. The modification of the Ainley-Mathieson method mostly 
affects the secondary and tip clearance loss coefficients, but has an 
affect on profile loss Yp for Mach numbers greater than 1, and a Rey- 
nolds number correction to the profile and secondary loss coefficients. 
3.2 The pressure loss coefficient Y is calculated in subroutine 

PL053. Figure (3-2) shows the sign convention used in this subroutine, 
where jangles are blade angles and OC angles are gas angles. The ang- 
les as shown are in the positive direction. The (') implies relative 
angles. This sign convention and notation applies only to subroutine 
PL03S in order to have a standard form for stator and rotor. 

Figure (3-3a) and 3-3b) coupled with equation (3.2-1) is used 
to calculate the profile loss coefficient Yp for 0 incidence, where 
$ 1 = 0 values come from figure (3-3a) and jQj = 0^Z comes from figure 
(3-3b) The pitch to chord ratio (P/C) is the argument for entering 
figures (3-3a) and(3-3b). The correction for thickness to chord 
ratio (t/C) is done in equation (3.2-1) along with the correction for 
the ratio of inlet blade angle to exit gas angle (/^i/^). The value 
obtained is Yp for zero incidence (i=*0). 

For off-incidence calculations figures (3-4) a,b,c, and figure 
(3-5) are used. Figure (3-4) a,b,c determine the stalling incidence 
so that figure (3-5) can then be used to correct Yp for of f- i ncidence 



cases . 
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Yp - j^p(/5,-0) + (£//06) 2 [*p(<3- OG)-YpCS/-oDJ (3.2-1) 
ft/c] (A/OQ 

L9* 2 J 



O^P/C=.75) 



* 2 . 

(OC2/(64 z (P/C-.7 5)) 



(3.2-2) 



Figure (3-4b) is entered with the argument of pitch to chord 
ratio (P/C) • The ratio of exit ?as angle (X^to exit angle£/£P/C=. 75) 
is then determined. It is . now possible l it o determine 75) Having 

de termi ned(y^ P/ C*. 7 5) , fiigurs ; (3^6c>) IB aaw entered wi th 6^( P/C 1 *^ 
to get the stalling incidence i s (P/C«. 75). Figure (3?-4a) is now entered 
•/ ith exit gas angle d P/C to get ^ i^ = i^ - i^(P/C=.75). can 

then be' solved for. Having incidence angle i and i^ the ratia 
(i/i^) for entering figure (3-5) to determine ratio Yp/ Yp( i»0) ds known. 
Having previously solved for Yp(i«0) in equation (3.2-1) Yp is now 
known for off incidence. 

Ys = .0334(C/H)(cosa£/cos^Z (3.2-3) 



For secondary losses Ys is computed by equation (3.2-3) 
from reference (8). Where Z is the Ainley loading parameter from 
reference (2). 



Z 




cos ^ pCz 
c°s 3 



(3.2-4) 



06" ” ( (y^ — / 2 (for stator) (3.2-5a) 

(tfn - (O^C fc/)/ 2 



(for rotor) 



(3. 2- 5b) 



-20 
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-1.2 -0.8 -0.4 0.0 0.4 0.8 1.2 

(c) 



Figure 3-4. Off-design perfornnance of cascades of turbine 
blades. Re - 2x10^- .(a)variation of stalling incidence 
with nitch-chord ratio;(b)variation of/l^ith p/c ; (c)s tal 1 - 
Inc Incidence of turbine blade sections when p/c - 0.75.^ 
From D. j.Ainley and G.C. R.Mathleson(2)) 
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^Wi 



Figure 3-5. Variation of relative profile loss with 
relative incidence. (Fron D.G.Ainley and G.C. R. 
Mathieson(2)) 




2(tan04 + can Of?) 



(3. 2 -6a) 



[?7^) “ 2 < ta "04+ tan <V ( 3 . 2- 6b 

For tip clearance losses fk is computed by equation (3.2- ) 
trora reference (8). 



Yk - il( C/ II) ( k/C) * L 



(3.2-7) 







Figure 3-0. effect of trailing edge thickness 
on oladc coefficients. ( From D. 3.Ai nley and 
G.C.ft.mathieson(2)) 



For equation (3.2-7) 

3 - .37 for shrouded blades 
3 - .47 for plain tip clearance 
3 « .00 for stator blades 
k - tip clearance 
C * chord 
11 - blade height 

having determined and fs they can be corrected for .Mach 
number and Reynold’s number from equation (3.2-3) and (3.2-9) 
respect i vc 1 y. 











% 
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(Yp) corrected ® Yp( l+60(Mn- 1 ) £ ) 
for Mn 

(Yp + Ys) corrected = ( Yp + Ys) 



ReJ "' 2 

1 



2x1 



(3.2-8) 



(3.2-9) 



Now the total loss coefficient Yt can be computed. by equation 
(3.2-10). This coefficient is now corrected for thickness of trailing 
edge to pitch ratio (te/P) from figure (3-6) to get the total loss 
coefficient for the blade row. Illustrations of an example calcula- 
tion can be found in reference (2) and (12). 



Yt - Yp + Ys + Yk 



(3.2-10) 
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3.3 In this section it will be explained how the loss coefficient 

Yt solved for in section 3.2 is coupled with compressible flow tables 
to estimate performance. The method explained is that outlined in 
reference (12). Another method is explained in reference (2). Cal- 
culations will be carried out only at the mean radius which should 
give a fairly accurate assessment of .the turbine performance. 

The calculations are started at the inlet of the turbine and 
worked through each stage where the exit conditions at one blade 
row is the inlet to the next blade row until the turbine exit of the 
last stage is reached. 

For inlet conditions the values for Toi and Poi will remain 
constant, the mass flow and corrected speed (n/ yjToi)/ ( N /jloi )dp will 
will be varied. This will allow the values of efficiency verses 
pressure ratio and 
culated as illustrated in Figure (3-7). 

3.3a Flow through the nozzle. 

(1) Select the mass flow. 

(2) Assume an exit Mach number. 

(3) With known blade angles subroutine S TOC HA is entered as 
in section (2.2) only this time the blade angles will be used to det- 
ermine the exit gas angles. 

(4) From compressible flow tables for the assumed Mach num- 
ber the values of (t^/Poz) and ( / N |To2j are determined, which will 
then determine V ^ anc * ^2 



> Will dll W w L 

n4Toi\ / (n [Tob i 

Poi )/ y Poi /c 



verses pressure ratio to be cal- 
'dp 



^2 (V 2 /{T52) 



(3.3-1) 
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IssentropU: EXXiclcucx vs Pressure ratio 

N/^ToV- 






(a) 



Po i / Poe 



7L» 

(b) 




Figure 3-7. Performance Curves for Axial Turbine (a) 
Issentropic efficiency vs pressure ratio, (b) Static 
efficiency vs pressure ratio, (c) Corrected nnass flow 
vs pressure ratio. 



Static efficiency vs Pressure ratio 
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/-I 

T 2 * T02XP 2 /P°2.) (3.3-2) 

(5) Subroutine PLOSS is now entered to get Yt (pressure 
loss coe ffiecient). 

(6) P 02 . can be determined from known values of Yt and 
(P 2 /Po2. ) from equation (3.3-3) 

Po^ = Pol/(Y t (l-P 2 /Poz.) + 1) (3.3-3) 

(7) The mass flow parameter Q2 can now be determined. 

Q2 = ,^ To (3.3-4) 

A 2 Po 

where A 2 = An 2 Cos£X> (3.3-5) 

(8) From the compressible flow tables for the value of Q 2 

the Mach number can be determined, and if this value is close to the 
assumed Mach number the inlet conditions to the rotor can be solved 
for the new values of V Mach number, and T 2 will be used 

to iterate the process until the assumed Mach number approximately 
equals the calculated value. Also as the Mach number approaches (1) 
a calculation will be performed to see if the row is choked. (ie; 
maximum flow for given parameters). This will be covered in section 
3.3c. 

Now that the exit Mach number has been determined the 
exit conditions of the nozzle can be determined, which will be the 
inlet conditions to the rotor. Where is determined from equa- 
tion (3.3-1)^ T 2 from equation (3.3-2) and ?2 from compressible 
flow data. 



P 2 - Poz(P 2 /Poz.) 



(3.3-6) 
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3.3b Flow through the rotor 

(1) WithC^and V ^known the relative rotor angle can 



be calculated. 



tan (iS 2 = tan^ 2 - 



Umo 



where Um = % 



V2COSA2 



(3.3-7) 



and determine the incidence angle (i,) where is the blade angle 

and i - ^2 



(2) The relative velocity W 2 and relative total temper- 



ature, can then be calculated. 



, 2 9 

W 2 - ( V 2 cos [X2 ) ( H* tarr ^ 2 ) 

2 

To 7-r e 1 ““ T> + W 2- 



2go JCp 



(3.3-8) 

(3.3-9) 



(3) W 2 /^To2. re ^ can be calculated, which will be used to 
enter the gas tables to get P 2 /PoZ re i which will allow Po2. ^ 

to be calculated. 

(4) The same procedure of assuming exit Mach number and 
reiterating until the assumed Mach number equals the calculated one 

is followed. P°3 re p <^3 > W^t^rel ,and r^rel'Will then be known. From 
this relative Mach number (P 3 /Po 3 re i) can be determined, which will 
enable P^ to be solved for. The absolute parameters OCy* 

P 03 , T^ and T 03 can then be calculated. 



tan^l 3 = tar ^ 3 + 
where Um » r^CO 



Um 



W 3 co ^ 3 



(3.3-10) 



I 

j 
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(W 3 cos£3) (1 + 



tano^-3) 



(3.3-11) 



r°3rel 



W2 

2<>oJCp 



/ 

P 03 = P3(To 3/T 3 ) yTT 



(3.3-12) 



(3.3-13) 



(5) If this is the last stage, e f f i ciency can be calculated 
by equations (3.3-14) and (3.3-15). If not the exit conditions will 
determine the inlet conditions to the nozzle of the next stage. This 
process will be iterated until the parameters have been determined for 
every stage in the turbine. 
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tt 



(Toi 

Toi 



Toe) 



1 - (Poe/Poi) 



rr-1 
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ts 



(Toi 
Toi l 



Toe) 

£l 

(Pe/Poi) * 



J 



(3.3-14) 



(3.3-15) 



3.3c Determination of whether a row of turbine blades have choked 
is carried out when the Mach number approaches (1). For the computer 
program this process will be carried out when the Mach number is great- 
er than 0.9. 

(1) Select a slightly higher value of Q ( HLl£2 — ) calling it 

APo 

Q* • 

(2) With this Q* enter the compressible flow tables to get 
corresponding values of P/ Po and solve for Poj n /Po 0ut assuming Yt 



rema i ns constant . 
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(3) Now calculate the new mass flow rn 1 . If this is less than m 
the flow has choked, Yt and Q wilL be assumed to remain constant with 
the flow going supersonic at the exit of the blade row. If in' is not 
less than m calculations will be resumed as was illustrated in section 
(3,3a) and (3.3b). 



From equation (3.3-3) determine Pout. With equation (3.3-4) Aout can 
be determined. Since Aout = ^annulus cos ^ be solved for in 
equation (3.3-16). The calculations are then resumed in the same 
manner as in section (3.3a) and (3.3b). The validity of this method 
depends on whether Q and Yt remain constant after the turbine blade 
row has cnoked. 



(4) Enter compressible flow tables with Q for supersonic flow 



to get corresponding values of ( P 0 ut/ Po out ) » 





(3.3-16) 
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Chapter 4 

Operation of Computer Program 

4.1 Chapters 1, 2, and 3 described how the calculations are 
performed in the computer program. This chapter will show how the 
flow of information is passed through the computer program to get 
the desired results. This information is not needed to use the 
computer program since that is covered in Appendix (2). 

4.2 The steps of calculations in the computer are in the same 
sequential order as illustrated in chapters 1 through 3. This is 
shown in the simplified block diagram figure (4-1) 

4.3 The method used in the computer program to interpolate 
curves is a 4 point Lagrange Polynomials with a worse case of a 
linear interpolation due to break down of Lagrange Polynomials for 
monatomic increasing or decreasing function. Subroutine FIG and 
BK are used to do this interpolating. Subroutine FIG selects the 
4 data points closest to the argument X for passing to subroutine 
3K which actually does the interpolating to get the ordinate Y as 
an output. If the argument exceeds the limits of the data points 
the value closest to the argument X is used to get the ordinate Y. 
The reason this method was chosen over finding formulas by least 
squares or simular methods is this method can be adapted to chan- 
ges in data without solving for new formulas or coefficients. The 
Lagrange Polynomials uses a monatomicly increasing base to get the 
ordinate Y. If the base is not monatomicly increasing such as 



SIMPLIFIED BLOCK DIAGRAM OF COMPUTER PROGRAM 
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values of Mach Number vs Ps/Po, subroutine FIG rearranges the base so 
that the base is monatomiciy increasing and the ordinate values are 
in corresponding relationship with the base. For values of base such 
as (m^Toi )/Apoi ) which both increase and decrease the data ooints are 
divided into two parts with the part that increases as one set of data 
points and the part that decreases as another. The ordinate does not 
have to be a monatomic function. The problem of increasing and 

decreasing bases can be encountered in the performance estimations be- 
cause the oase and ordinate are i nterehan^ed in some of the steps. 
Figure 4-2 illustrates how subroutine FIG and 3K are used to inter- 
polate data. 

4.4 Figure (4-3) illustrates the method the computer program uses 
to determine the flow area at station 1. In most cases for a value of 
epsilon (6) = .001 the convergence is reached in less than 10 iter- 
ations. This method is used for solving for the exit flow area and 
also for station 2. The method for solving the other parameters 

were will documented in chapers 1 thru 3 and will not be repeated 
here. 

4.5 From output generated by the computer program, it is felt that 
the computer program can be very useful from early preliminary design 
to start of detailed design, since it has the versatility to accept 
more input data as the designer firmly sets the characteristics of the 
Turbine. As for cooling requirements it is felt that compensations 

in the input data could account for losses attributed to cooling. 



-K 



BLOCK DIAGRAM OF SUBROUTINE FIG AND 3K FOR INTERPOLATING 
FUNCTIONS £(x) from MAIN PROGRAM OR OTHER SUBROUTINES 
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Figure 4-2. Flow / Mow, 3-52 % j 
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APPENDIX I 



Symbols and Notation used in this Text. 



SYMBOLS USED IN TEXT 



A Area 

C Chord of turbine blade 

Cp Specific heat (btu/lbm°R) 

go Constant in Newtons Law 32.174 ( f tl bm)/ ( 1 bfse 

h Blade height 

ho Stagnation Enthalpy 

H/C Aspect ratio 

i Incidence angle 

I Moment of inertia 

J Energy conversion factor 778.16 ( ftlbf )/ (btu) 

k radius of gyration 

M Mach Number 

M crit Vx/( ^7f 50R T 0) k 
m Mass Flow 

N Number of Stages 

n Number of Stage 

Po Stagnation Pressure 

P Static Pressure 

P/C Pitch to Chord ratio 

R Gas cons tant(ftlbf )/(lbm £> R) $^JCp 

Q Corrected mass flow m^To/PoA 

To Stagnation Temperature 

T Static Temperature 



To 



Thickness of Outer Rim 
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Ti 

SM 

U 

V 

Vx 

W 

Wo 

Wi 

Y 
R 
r 
t 

t/C 

Z 

Zi 

Zo 

Z(r) 



cx 




7L 



Thickness of Inner Rim 

Non Dimensional Section Modulus 

31ade Speed 

Velocity 

Axial Velocity 

Relative Velocity 

Width of Outer Disc Rim 

Width of Inner Disc Rim 

Pressure Loss 

Radius 

Radius 

Thickness of Turbine Blade 
Thickness to Chord Ratio 
Section Modulus 

Thickness of Mid Section of Turbine Disc at 
Thickness of Mid Section of Turbine Disc at 
Thickness of Disc as a Function of Radius 
Absolute Gas Angles (except section 3.2) 
Relative Gas Angles (except section 3.2) 
Stagger Angle 

Ratio of Specific Heats Gamma 
Ef ficiciency 
Stress 
Dens i ty 

Psi Loading Coefficient 



Inner Radius 
Outer Radius 



67 



J2 


Reaction 
Speed (rpm) 


CO 


Speed (rad/sec) 



Subscripts 



i 


Inlet 


e 


Exit 


t 


Total 


s 


Stage 


P 


Poly tropic 


s 


Static 


t 


Tip 


m 


Mean 


h 


Hub 



c 



constant 
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APPENDIX II 



User's Guide To the Computer Program with Comparison of 
Turbine Designed with this Computer Program to 
one which was.. Designed by Airesearch. Sample 
Input and Output of Turbine 
Designed by Computer 
Program 
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II-l The computer program can be used to calculate physical and 
thermodynamic properties of an axial flow turbine.. It has an optional 
ability to do performance calculations. The optional phase is very 
expensive to operate since it does both design and off design calcu- 
lations. From operating the computer program and comparing the cycle 
calculations with performance calculations the cycle calculations 
are within 37. of performance calculations at oesign point. For the 
early preliminary phases of turbine design it is recommended that 
just the cycle calculations and mechanical design phases be carried 
out with the computer program punching output cards to have data 
to perform the performance calculations at a later time 

A:.second program which reads in the punched output cards in the 
same order as they were punched in the first program is provided to do 
only performance calculations. This allows the performance estimations 
to be carried out with runniding the complete program again. 

For future references program #1 will refer to the complete 
computer program while program #2 will refer to the program which only 
does performance estimations. The complete computer program can not 
be operated in computers which have less than 210 k of memory. 

II-2 Pages 71 to 76 show how the computer cards are to be punched. 

The 1st card is a label which will be printed out in the same content 
as it was entered on the computer card. Cards 2 thru 4 are the input 

data cards of known turbine parameters. Cards 5 thru NNN are the 

compressible flow tables and viscosity of the gases with corresponding 
temperatures. Starting from card 5, 1st the Mach Numbers are read in 
, 2nd the corrected mass flow m |to/ PoA, 3rd the pressure ratio P/Po^ 
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4th the corrected velocity V/^To, 5th the temperatures which correspond 
to the viscosity data, 6th the viscosity data. Page 77 is a set of 
compressible flow tables and viscosity data for air. 

II-3 Pages 73 thru 32 are comparison data of computer program designed 
turbine with one which was designed by Airesearch reference (20). As can 
be seen from the data some of the parameters are different because of 
lack of knowledge of some of the design parameters which was used in the 
design of the Airesearch Turbine, also some simplifications are 
incorporated into the computer program for solutions of output para- 
meters. 

LI-4 A list of sample output and input is provided from page 83 to 
93 for turbine which was compared with Airesearch turbine. As can 
be seen on page 91 of output data, the performance data is very sim- 
ilar to the cycle calculations at design point. This justifies the 
argument of doing only cycle calculations and mechanical calculations 
when off design data is not required. Figure II- 3 is a plot of the 
data points from performance calculations. Where some of the curves 
are discontinuous is caused by a limit of 10 iterations for solution 



at a blade row. . 



Input Data Cards 
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* (Must be specified) 

# (If value of parameter is not known place value to right of # 

on computer card) 



Card 1 

Column 1-80 



Card 2 

Column 1* 



Column 2 * 



Column 3 and A 
Column 3 * 



Alphanumeric label (ie, M 3as Tur- 
bine number l 11 ) 

# Blank Card 

Ml (Integer) 

0 compressible flow data to be 
entered. 

1 No compressible flow data to be 
entered. 

M2 (Integer), 

0 Power in horse power specified. 

1 Pressure ratio specified. 

Blank 

M (Integer) Type of turbine. 

1 Constant hub 

2 Constant mean 

3 Constant tip 



Note Where instructions inform user to place dummy values on 
cards these values are used to inform computer that they 
are not specified and to perform optional calculations to 
determine the parameter. 



72 



Column 6-10 


KLM (Integer right justified) 


(leave blank i f Ml s 


* 1) Number of compressible flow el- 

ements to be entered. (ie num- 
number of Mach Numbers to be en- 


Column 11-15 


entered) ( Maximum 100) 

KIM (Integer right justified) 


(leave blank i f Ml = 


= 1) Number of Temperatures at which 

values of (viscosity) x 10 ^lbm 

seC ft. 

is defined. 


Column 16-20 * 


Gamma (F5.3) 


Column 21-25 * 


CP ( 3tu/lbm °R) 
(F5.3) 


Column 26-35 


Power (Horse Power) ( F10.1) if 


(lbm/sec ft) 


M2 ■ 0 this must be specified 
# 9000. 


Column 36-45 * 


W (mass flow) (lbm/sec) ( F 10.3), 


Column 46-50 * 


TC (thickness to chord ratio) (F5.31, 
normally about # 0.2 


Column 51-53 * 


TET (thickness of trailing edge to 
pitch ratio) (F5.3) normally about 
# 0.02 


Column 56-60 * 


TIPCLA (Tip clearance) (in) (F5.3) 
normally about # 0.015 



73 



Card 3 




Column 1-5 * 


p PSI (F5.3) 

Normally between 1.5 and 2.9 
# 2.0 


Column 6-10 * 


RTE (Max Tip raduis in ft at exit) 
(F5.3) 


Column 11-15 * 


Ratio T (Rte/^he) Ratio of Tip to 


Column 16-20 


hub at exit (F5.3) # 1.5 

(J\hub (exit) Reaction at hub (F5.3) 

normally between .00 to 0.4 # .05 


Column 21-30 * 


OMEGA (RPM) (F 10. 0) 


Column 31-32 


KLMl (INTEGER) Right justified if 
Compressible flow data provided 
the value at which 

TbcA 

is a maximum in (for 

sample problem KLM 1 « 50) (Maximum 70) 


Column 33-39 


Blank 


Column 40 * 


M4 (Integer) 

0 performance calculations to be per- 
formed. 

1 No Performance Calculations to be 
performed. 


Column 41-45 * 


SPECWT (Specific Weight of turbine 
disc and blades (lb/in^) (F5.3) 



/) normally about 0.3 



Column 46-50 * 
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Y MOD (Young's Modulus psi) (f 15.0) 
Normally about # 30000000. 


Column 61-70 * 


Cl (Value of allowable blade stress 
in psi) (F 10.0) Normally about 
# 50000. 


Column 71-80 * 


C2 (Value of allowable blade steady 
stress in PSI) (F 10. 0) Normally 
about # 100000. 



Card 4 



Column 1-5 * 


ENT (Value of total efficiency 
desired. # .9 


Column 6-15 


Radius of Shaft in ft (F 10.5) # - 0. 1 
if not known 


Column 16-20 


P ratio (Pressure Ratio) (F5.3) 
i f M2 ■ 1 this must be specified 77 1.5 


Column 21-30 * 


ULTSRE (Ultimate Tensile Strength of 

Material of Disc (psi) Value for disc 

stress is calculated from this. (F10.0) 

avg = .75 ULTSRE 
2 

(1.2) 

Normally about # 130000. 


Column 31-35 


ALPHAH (Exit angle in radians) (F5.3) 
Normally between (0.0 to - .3) # - .1 


Column 36-45 


EMACH - Vx Normally between 

soaTo 

.2-1.0 #-.5 



Column 46 


M3 (Integer) ^5 

0 Normal Blades 




1 Shrouded Rotor Blades. #o 



Column 47-48 * 





NKL (Number of turbines to be de- 
signed by computer program. (In- 
teger right justified) All sets of 
cards except for last turbine des- 
ign will have 0 in column 48. The 
last set of cards will have the nGm- 
ber of sets that were entered. This 
allows the computer to stop after com- 
pleting the last turbine design. 


Column 49-50 


N (Number of stage Integer right jus- 
tified) (Maximum 10) If this value 
is not known, enter 0 and computer 
will solve for number of stages. #0 


Column 51-55 * 


AREARA ( F5.4) 

Tip to base area ratio of turbine 
blade normally between .25 to .35 

# .25 


Column 56-60 


TOO (Gives radial thickness of outer 
rim of disc. by To-** Wo/TOO normally 



between 2 to 4 (F5.3) # - 5. 



Cards 5-NNN 
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These cards are optional. If performance calculations are 



to be performed they must be furnished. If more than one turbine 
is to be designed compressible flow data must only be furnished 
for 1 st turbine and the value of Ml on card 2 for each following 
turbines will be "l". If the compressible flow parameters are dif- 
ferent a set must be furnished. A set of compressible flow data and 
viscosity data for air are provided with the computer program. 

Card 5 - N RMAC (Mach Numbers) (20F4.2) 



same number as KLM 



Card (N + 1)-K 




PoA sec lbf 
Same numbers as KLM 



Card (K+l ) -J 



PSPT (Ps/Po) (10F8.5) same 



number as KLM 



Card (J+l)-L 



VELTOT (V/J To) (10F8.5) (ft/sec) 

T°o 



V « 



same number as KLM 



Card (L + L) - LL 



TEM (°?J (20F5.0) same number as 



KIM temperature values that cor- 



respond with viscosity data 



Card (LL + 1)-NNN 



V (Viscosity of fluid ( 1 bm. sec/ f t ),x 
10 7 ) (20F5.0 ) 
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dcmprt^sibft^Bou/TaUts for Cr>n)b' 

-•' tbr» \°R •' G+/$eZ)Aon 

777 It,*- ' K 



MACH NO 



•5«C lb? 

W*SQRTJTI/(A*P) 



Tab/ti f»y r nn ,b^-ftO" P rod ucts, 

nz-feSQuce C 2 -) 



_ 0. 02000 
0.04000 
0.06000 
0.08C00 
0.10000 
0.12000 
__0. 14000 
0.16000 
0.17000 
0.20000 
0.22000 
0.24000 
_ .0.26000 
0.28000 
0.30000 
0.32000 
0.34C0O 
0.36000 
0.38000 
0.40000 
0.42000 
0.44000 
0.46000 
0.48000 
_0. 50000 
0.52000 
0.54000 
0.56000 
0.58000 
0.60000 
0.62000 
0.64000 
0.66000 
0.68000 
0.70000 
0.72000 
__0. 74000 
0.76000 
0.78000 
0.80000 
0.82000 
0.84000 
0.86000 
0.88000 
0.90000 
0.92000 
0.94000 

0. 96000 
_ 0.98000 

l. 00000 
1.02000 
1.04000 
1 .06000 
1.08000 

1 . 10000 
1.12000 
1.14000 
1.16000 
1.18000 
1.20000 



0.01811 
0.03J96 
0.05388 
0.C7171 
0.08945 
0.10706 
0.12448 
0.14154 
0.15885 
0.17562 
0. 19226 
0.20849 
0.22446 
0. 24049 
0.25612 
0.27135 
0.28590 
0.30026 
0.31435 
0.32763 
0.34118 
0.35386 
0.36627 
0.37821 



PS/PT V/SQRTm 

0.99973 0.96 15 

0.99893 1.91184 



0.99760 
0.99574 
0.99335 
0.99040 
0.98700 
0.98317 
0.97870 
0.97380 
0.96835 
0.96260 
0.95630 _ 
0.94950 
0.94200 
0.93440 
0.92650 
0.91800 
0.90920 
0.90C30 
0.89080 
0.88100 
0.87060 
0.86000 



2.86589 

3.82740 

4.77028 

5.72284 

6.67094 

7.61C08 

8.56787 

9.50329 

10.45287 

11.36668 

12.31328 

13.26734 

14.22884 

15.14564 

16.04750 

16.99411 

17.90344 

18.82024 



~£>r y-/333 



it 






Jb*A 



18.82024 j * 

19.75192 07 

Air- 'Jron %C-ftfettCC (l4) 



0.38948 „ 0.84950. 
0.40061 0.83860 



0.41108 

0.42114 

0.43067 

0.43979 

0.44838 

0.45669 

0.46421 

0.47159 

0.47796 

0.48433 

0.49010 

0.49560 

0.50030 

0.50459 

0.50815 

0.51170 

0.51479 

0.51734 

0.51948 

0.52123 

0.52257 

0.52337 

0.52391 

0.52418 

0.52391 

0.52337 

0.52257 

0.52136 

0.51989 

0.51801 

0.51599 

0.51358 

0.51076 

0.50795 



0.82720 
0.81570 
0.80400 
0. 79220 
0.78C20 
0.76780 
0.75550 
0.74270 
0.73070 
0.71800 
0.70520 
0.69140 
0.67910 
0.66650 
0.65500 
0.64080 
0.62760 
0.61410 
0.60250 
0.58900 
0.57550 
0.56500 
0.55100 
0.54000 
0.52800 
0.51700 
0.50400 
0.49190 
0.48000 
0.46820 
0.45650 
0.44530 
0.433.30 
0.42250 



20.63144 
21.54077 
22.45010 
23.33708 
24.27625 
25. 16321 
26.05019 
2 6.95207 
27.80922 
2 8.65892 
29.54590 
30.41052 
31.29750 
32.09502 
32.94472 
33.80190 
34.65160 
35.47894 
36.29883 
37.05164 
37.93861 
38.76595 
39.60822 
40.39828 
41.14365 
42.03061 
42.70809 
43.41698 
44.16234 
44.98222 
45.69031 
46.36113 
47. 17357 
**7.85104 
48.59720 
49.34256 
50.02827 
50.74382 
51.39229 
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COMPARISON OF COMPUTER PROGRAM TURBINE DESIGN WITH 78 

AIRESEARCH 9000 HP FREE VORTEX AXIAL FLOW TURBINE 



Inlet Condi t ions 






Gas 


Ai r 




Poi (psi ) 


312.5 




Toi ( °R) 


1767.8 




m(lb/sec) 


119.421 


Design 


Parameters 






i> 


1.8 




(X® 


0.0 


Assumed 


Design Parameters 






0.9168 




Cp(btu/lb°R) 


.2715 




Gamma 


1.333 



Results 

Figure I I- I is a comparison of the turbine wheel dimensions. 



Figure II-2 is a comparison of design point performance calcu- 
lations of total and static efficiencies as a function of tip clearance 
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COMPARISON OF OUTPUT DATA FROM COMPUTE 
AIRES EARCH GAS TUR3INE 



Parameter 

Toe (°R) 
M cri %t exit 
Poe (psi) 

hub(statioa 1) 

radius(£t) 

& l (des) 
C^ l (de§) 

tip(station 1) 

radius(ft) 

/3 i (de - > 
(yj ( de s) 

Huo(station 3) 

radius( f t) 

& 3 ^ de s) 

( de o) 

Tip(station 3) 

radius( f t ) 

/3 3 ( d e §) 
OC 3 < de S> 
d cf(P si > 

Wheel Weight(lbf) 



Airesearch 

1571.4 

186.0 

Velocity Diagrams 

.4383 

-62.02 

0.0 

.530 

-67.15 

0.0 

.4383 

-60.95 

0.0 

.6215 

-68.26 

0.0 

37,400. 

18.7 



R PROGRAM WITH 

Computer Program 
1554.0 
.292 
177.3 

.429 

-63.96 

0.0 

.429 

-70.93 

0.0 

.429 

-63.96 

0.0 

.636 

-72.99 

0.0 

37,396. 

21.46 



30 



Conclusions : 

Most of the descrepancies in parameter values can be attributed 
to different input values for the computer model. This is caused by 
not knowing the assumptions or definitions that were used in designing 
the Airesearch Gas Turbine. From this limited comparison it is felt 
that the computer model can be used for designing marine gas turbines 
since it has the versatility to change parameters which are attributed 
to specific turbines. For high performance turbines (aircraft) it ' 
should be compared with known aircraft turbines to verify the validity 
of the computer program and adjust the variable parameters to match 



them. 
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Figure IL-2, Comparison of Efficiencies verses 
Tip Clearance of Performance Calculations. 
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with Computer Program, (a) Issentripic Efficiency (b) 
Static Ef f iciency. (c) Corrected mass flow. 
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APPENDIX III 

Computer Programs developed for the Design 
of a Free Vortex Axial Flow Turbine. 

Computer Program 1 Complete design 

Computer Program 2 Performance Calculations Only 
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Preliminary design of 
a free vortex axial 
flow turbine. 



